HYPERBOLIC BETA INTEGRALS 
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Abstract. Hyperbolic beta integrals are analogues of Euler's beta integral in which the 
role of Euler's gamma function is taken over by Ruijsenaars' hyperbolic gamma function. 
They may be viewed as (q, (f)-bibasic analogues of the beta integral in which the two bases 
q and q axe interrelated by modular inversion, and they entail g-analogues of the beta 
integral for \q\ = 1. The integrals under consideration are the hyperbolic analogues of the 
Ramanujan integral, the Askey- Wilson integral and the Nassrallah-Rahman integral. We 
show that the hyperbolic Nassrallah-Rahman integral is a formal limit case of Spiridonov's 
elliptic Nassrallah-Rahman integral. 



1. Introduction 
Euler's gamma function is defined by 



PCX) 

T(z) = / x z - l e~ x dx, Re(z) > 0. 
Jo 



In the fundamental paper JH], Ruijsenaars defined gamma functions of rational, trigono- 
metric, hyperbolic and elliptic type, Euler's gamma function being of rational type. Accord- 
ingly, one expects to have extensive theories on special functions of rational, trigonometric, 
hyperbolic and elliptic type. The rational and trigonometric special functions are precisely 
the special functions of hypergeometric and basic hyper geometric type, which have been 
thoroughly studied (see e.g. [2] and [IB])- Systematic studies of hyperbolic and elliptic 
special functions have commenced only recently, see e.g. [20] > [HI] for the hyperbolic 
case and [12], [23, [21] for the elliptic case. 

A basic step in the development of special functions of a given type is the derivation of 
the associated beta integrals. The beta integral of rational type is Euler's beta integral 

(1.1) I' x a -\l - xf- l dx = ffl^fj , Re(a),Re(6) > 0. 

Jo r ( a + b > 

Trigonometric and elliptic analogues of beta integrals (involving the trigonometric and the 
elliptic gamma function respectively) have been studied in detail. The goal of this paper 
is to derive hyperbolic analogues of beta integrals. 

The importance of beta integrals lies in its variety of applications. The rational beta 
integral is the normalization constant for the orthogonality measure of the Jacobi polyno- 
mials and, more generally, multivariate analogues of the beta integral arise as normalization 
constants in the theory of zonal spherical functions on compact symmetric spaces. Beta 
type integrals also appear in number theory, combinatorics, conformal field theories and in 
certain completely integrable systems. 

l 
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Trigonometric beta integrals have been studied extensively, see e.g. 0, [7], 
[To] . In this case the role of Euler's gamma function, as well as of functions like x a ~ l 
and (1 — are taken over by (quotients of) Ruijsenaars' ^Hj trigonometric gamma 

functions, or equivalently, (quotients of) g-Pochhammer symbols. The g-Pochhammer 
symbol is defined for |g| < 1 by 

oo 
j=0 

and products of g-Pochhammer symbols will be denoted by 



(a 1 ,...,a m ;q) oo = Y[(a j ;q) c 



;i. 2) 



in 



{ai4\---^mzi 1 ;q) oo = Y[(a j z j ;q) oQ (a J z j 1 ;q) oQ . 

3=1 

A trigonometric beta integral which closely resembles (jl.lj) is the contour integral 

, s If {~q* z~ l , -qh; q) ^ dz (qa,qb;q) c 

(l.dj 



2tt« Jt (-q^az- 1 , -qhz; q)^ z (g, qab; g)^ 
for |a|, |&| < \q~^\, where 

T = {z G C| \z\ = 1} 

is the positively oriented unit circle in the complex plane. We call (jl.3j) the Ramanujan 
integral, since it is closely related to one of Ramanujan's trigonometric beta integrals from 
his lost notebook (see jl] and references therein). A second trigonometric analogue of 
is the Askey- Wilson integral [TJ, 

(1 4) J_ f ( z±2 ^)oo dz = 2(t 1 W 4 ;g) 00 

for parameters tj G C with \tj\ < 1 (j = 1, . . . , 4). The Askey- Wilson integral is the normal- 
ization constant for the orthogonality measure of the celebrated Askey- Wilson polynomials 
[7j. The Nassrallah- Rahman jTH] integral is the trigonometric beta integral 

(L5) 1 f (z ±2 ,Az^q) oQ dz mU( At ^), 



2m JT n.tofe^* 1 ; ?)oo Z (?5 ?)oo n <*<m<4(*fc*m; 

with A = to^i^2^3^4 and with parameters tj G C satisfying \tj\ < 1 (j = 0, ...,4). The 
Nassrallah-Rahman integral is the normalization constant for the biorthogonality measure 
of a five parameter family of io0g rational functions, see ^H] • Note that the Askey- Wilson 
integral ()1.4|) is the special case to — of the Nassrallah-Rahman integral (|1.5|) . The 
Nassrallah-Rahman integral is the most general trigonometric analogue of Euler's beta 
integral known up to date. 
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One way to introduce Ruijsenaars' ^H] hyperbolic gamma function, which will take 
over the role of the g-Pochhammer symbol in the hyperbolic analogues of (jl.3|) . (J1.4)) and 
fll.f)|) . is by its explicit infinite product realization, see |22|, 120] • Explicitly, we define for 
r G C \ {0} two deformation parameters q = q T and q = q T by 

(1.6) q = exp(27ur), q = exp(— 2iri/r). 

The two deformation parameters are thus related by the transformation r i— > — i, which is 
one of the modular transformations preserving the upper half plane 

U= {z e C|Im(z) > 0}. 

We now define for r G H (so \q\, \q\ < 1) a (g, g)-bibasic analogue of the Pochhammer 
symbol called the r-shifted factorial by 

(exp(— 2iriz): o) 

Shintani's [22] crucial result implies that [z; r]oo extends continuously to r G M <0 for generic 
z G C (in which case \q\ = 1 = \q\). The resulting function, which we still denote by [z\ r]^, 
depends meromorphically on z G C and can be expressed explicitly in terms of Barnes' 
double gamma function, or equivalently in terms of Ruijsenaars' jTH] hyperbolic gamma 
function or Kurokawa's double sine function (see [201 Appendix A]). In the appendix, §6, 
the interrelation with Ruijsenaars' hyperbolic gamma function is made explicit and relevant 
properties of the hyperbolic gamma function are translated to the r-shifted factorial [z; t]^. 

The hyperbolic analogues of ()1.3|) . ()1.4|) and (jl.5|) can now be formulated explicitly as 
follows. We define the shorthand notations 

m 

3=1 
rn 

[ai ±Zi,...,a m ± z m ; r]oo = Y\\ a j + z i'i r ]°°K' ~~ z f r ]°° 

3=1 

for products of r-shifted factorials. We use the convention that q u = exp(27riru) and 
q u = exp(— 2'Kiu/r) for u G C. The hyperbolic analogue of the Ramanujan integral ()1.3|) 
becomes 

g-Ag A [± + a, ^ + /3; t]oo 
[i + a + ^jrjoo 

for Re(ra), Re(r/3) < and Re (a — | + , Re(/3 — | + < 0, where is the branch 
of the square root on C \ M, <0 which takes positive values on IR> (this branch of the 
square root will be fixed throughout the paper). Here the deformation parameter r may 
be either from the interior of the second quadrant of the complex plane, so Re(r) < and 
Im(r) > (in which case |g|, |g| < 1), or from IR<o (in which case |g| = 1 = |g|). The rather 



9 ~l~ 2T Zj ">-r ZlT 



2t 



a + ±- + a + z± + ±+P-z;T 



dz 
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mysterious looking q and q powers in the right hand side of (jl.8|) arise from an application 
of the modularity of the Dedekind eta function, 



fe^oo q ^q™ 



t e H. 



-IT 



The proof of (jl.8j) will given in §3. With the same conditions on the deformation parameter 
r, the hyperbolic analogue of the Askey- Wilson integral (jl.4)l becomes 



1 ~ 1 



. ( 9) . [1 ± 2z; r]oo ^ 2<j -<w) ^ [r, + To + r,. + t,, - 3: r| 



OO 



H 4 j=1 [r j ±z;T] o0 v 73 ^ 1 lli<fc< m <4[ r fc + r ™ ~ ^ r ]°° 

for parameters Tj G C satisfying Re(rj) < 1 (j = 1, . . . , 4) and Re((ri+r 2 +T3+T4— 3)r) < 1. 
The proof will be given in §4. The hyperbolic analogue of the Nassrallah-Rahman integral 
()1.5|) with deformation parameter r satisfying the same conditions as for the hyperbolic 
Ramanujan integral 1)1.80 . reads 

(11Q) r°° [l±2z,q-4±^;r] 00 ^_ 2q~^q^ Y\) =0 [ a ~ r 3 ~ 



Yl j=0 [Tj±Z;T] oo \T^VT n <fc<m<4[ r fc + T m - ^ T ]°c 

where a = ro + ri + r 2 + r3 + T4, with the parameters Tj G C satisfy Re(rj) < 1 (j = 0, . . . ,4) 
and Re (a — t _1 ) > 4. The proof will be given in §5. The hyperbolic Askey- Wilson integral 
(11.9)) is formally the limit case Ihi(to) — > — 00 of the hyperbolic Nassrallah-Rahman integral 
ffl~TU|) . see Remark EH) 

The proofs are based on the observation that a hyperbolic beta integral is essentially 
the fusion of a trigonometric sum and a trigonometric integral. The technique of fusing 
trigonometric sums and integrals will be discussed in §2. The hyperbolic Ramanujan 
integral (J1.8j) is then derived from fusing Ramanujan's 1^1 summation formula with the 
trigonometric Ramanujan integral ()1.3|) . while the hyperbolic Askey- Wilson integral ([1.90 
is derived from fusing Bailey's summation formula with the trigonometric Askey- 
Wilson integral (jl.40 . The hyperbolic analogue of the Nassrallah-Rahman integral is more 
delicate: it is derived from fusing the trigonometric Nassrallah-Rahman integral with a 
formula expressing the bilateral very-well-poised as a sum of two very-well-poised 807 
series. During the fusion procedure the sum of the 807 series is corrected in such a way 
that it becomes summable by Bailey's summation formula, leading eventually to ()1.1UJ1 . 

The fusion technique also reveals a close connection between trigonometric beta integrals 
and Macdonald-Mehta type integrals. Roughly speaking, Macdonald-Mehta integrals can 
be obtained by fusing trigonometric beta integrals with the inversion formula for the Jacobi 
theta function (regarded as a bilateral sum identity). We will encounter several one- variable 
Macdonald-Mehta type integrals in this way. In particular, we derive in §4 Cherednik's 
jH] one-variable Macdonald-Mehta integral as a consequence of the Askey- Wilson integral 
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For elliptic beta integrals, the role of Euler's gamma function is taken over by Ruijse- 
naars' ^H] elliptic gamma function 



;i-H) T(z; Pll p 2 ) = H — 



Pi F2 

kr^ra 



where pi,p 2 £ C are two arbitrary complex numbers with \pi\, \p 2 \ < 1. Spiridonov 
proved the following elliptic analogue of the Nassrallah-Rahman integral, 

/ 112 v 1 f IlJ=0 F(t 1 Z ±1 ;p 1 ,p 2 ) dz _ 2 U <k<m<i T ( t ktm;Pl,P2) 



2m 



where A = tot^tst^, with the five parameters £j G C satisfying \tj\ < 1 and \pip 2 \ < \A\. 
Here we used the notations 



r(ai,...,a m ;px,p2} = JJr(a i ;pi,p 2 ), 

i=i 
f;t 

r(ai2f 1 ,. . . ,a m ^ 1 ;p 1 ,p 2 ) = JJr(a i ^-, a J ^~ 1 ;p 1) p 2 ) 

i=i 

for products of elliptic gamma functions. The trigonometric Nassrallah-Rahman integral 
is the special case of the elliptic Nassrallah-Rahman integral (jl.12)) when one of the 
deformation parameters pj is zero. We show in §5.4 that the hyperbolic Nassrallah-Rahman 
integral (jl.lOj) can also be formally obtained as limit of the elliptic Nassrallah-Rahman 
integral ()1.12j) . It is based on a remarkable limit from the elliptic gamma function to 
the hyperbolic gamma function due to Ruijsenaars [T^]. The rigidity of the elliptic theory 
seems to prevent the existence of elliptic degenerations, and consequently there do not seem 
to exist elliptic analogues of the Ramanujan integral and of the Askey- Wilson integral. 

Remark 1.1. After the appearance of this paper as a preprint, Erik Koelink kindly pointed 
out to me that Ruijsenaars recently obtained another proof of the hyperbolic Askey- Wilson 
integral (|1.9|) as a spin-off of his studies of Hilbert space transforms associated to relativis- 
tic hypergeometric functions in his recent paper [21]. Subsequently van Diejen and Spiri- 
donov derived multidimensional analogues of the hyperbolic Askey- Wilson and Nassrallah- 
Rahman integrals (jl.9|) and ()1.10|) by completely different methods in [IT]. Their mul- 
tidimensional hyperbolic Askey- Wilson integral fTJ Thm. 5] can be rederived with the 
methods of this paper by fusing Gustafson's Jl] multidimensional Askey- Wilson integral 
with van Diejen's JU] multidimensional generalization of Bailey's very-well-poised 
summation formula. 



Notations: For \q\ < 1 we denote 



{zq a -q) c 
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This reduces to a finite product when a G Z. Products of [z; qj will be denoted in the 
same way as for [z; g) , see fjl.2|) . The r+ i0 r basic hypergeometric series with base q is 
defined by 



a 1; a 2 , . . . , a r+ i 
bi,...,b r 



q,z 



E 

m=0 



ai, • . • , a r +i; 



(g,6i, ...,6 r ; q) 



rn m 
6 ■> 



\z\ < 1. 



The special case of a very- well-poised r +i4>r series is defined as 

r+lWr ( a li a 4j • • • > a r+li 9) 2 ) = r+l<Ar 



l l 

a 1; ga^ , — qa\ , 04, . . . , a r+ i 
1 1 



E^ 



a^, — , qai/a^, . . . , gai/a r+ i 
1 - aig 2m (01,04, ...,a r +i;?), 



_-ai (g,gai/a 4 ,...,gai/a r+ i;g) ) 
The bilateral basic hypergeometric series r ^ r with base q is defined by 



rlpr 



a±, . . . , a r 
bi,...,b r 



q,z 



E 



m=— 00 



ai, 



a r ; q) 



(bi,...,b r ;q) 



m m 
Z « 



bi ■ ■ ■ b r 



CL\ • • • CL r 



< \z\ < 1. 



Acknowledgments: The author is supported by the Royal Netherlands Academy of Arts 
and Sciences (KNAW). 

2. Folding and fusion of integrals. 

In this section we explain the principle of folding and fusion of integrals, which is fun- 
damental for the study of hyperbolic beta integrals. 

We take r e H, so that \q\, \q\ < 1 for the associated deformation parameters q and q 
(see Q1.6]) ). The (rescaled) Gaussian 

(2.1) G T (z) = q$ 

is analytic, zero-free and satisfies the quasi-periodicity conditions 

G T (z + 1) = q^q z G T (z), G T (z - r _1 ) = q~^ exp(-2niz)G T (z). 

Suppose that <f>(z) is a ^-periodic, meromorphic function and consider the associated func- 
tion 

4>{z) = <P{z)G T {z) = <P{z)q^ . 
By the ^-periodicity of 4>, we still have the quasi-periodicity 

(2.2) 4>{z - t~ 1 ) = q-^ exp(-27riz)(j)(z). 

We now force <j>{z) to become one-periodic by considering the series 

00 

(2.3) <P + (z)= + 
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(this resembles the construction of automorphic forms as Poincare series). We assume that 
the function </> is such that the series (12. 3|) converges to a meromorphic function. In this 
situation, the resulting bilateral sum <fi + {z) is thus an one-periodic meromorphic function 
satisfying 

(2.4) (p + {z-T~ l ) = exp(-2niz)(j) + (z). 

On the other hand, the Jacobi theta function 



(2.5) $_i(z) = G_ i (m) exp(2nimz) = q 2 exp(2Tcimz) 

m=~oo m=—oo 

is a one-periodic entire function which satisfies the same quasi-periodicity (|2.4|) . hence 

00 

(2.6) <p + (z) = i{z + m) = ®{z)-&_x{z) 



m=— 00 



with $ an elliptic function with respect to the periods 1 and — . 

If the bilateral series <p + {z) is entire, then the elliptic function is automatically a 
constant. This well known fact follows from the observation that the Fourier coefficients in 
the Fourier expansion 4> + {z) = J2 n a n exp(27unz) satisfy the first order recurrence relation 
a n+ i = q2 +n a n (n G Z) due to the quasi-periodicity (|2.4j) of </> + . 

Choose now in addition a one-periodic meromorphic function ip{z). In the examples 
treated in this paper, the integral 

(2.7) / i)(x)$(x)$_i(x)dx 

Jo 

over the period cycle [0, 1] will be closely related to the trigonometric beta integral (jl.3|) . 
(Oil or (JT3J). The integral 



(2.8) / ip(x)(j)(x)q 2 dx = / ip(x)(j)(x)dx 

J —00 J —00 

may then be viewed as the fusion of the integral (|2.7J) and the sum ([2.6)1 . since folding the 
integral and using the one-periodicity of ip implies 



/OO 2 f 1 

il)(x)(j)(x)q~ dx = / 
-00 „ Jo 



i/j(x)4>(x + m) (ix 



m=— 00 
•1 



ijj(x)4> + (x)dx = / ■?/'(x) ( l ) (x)'(9_ 1 (x)dx 



(2.9) 



by the bilateral sum (|2.6p . In all applications below, (j2.9J) can be justified by a straightfor- 
ward application of Fubini's theorem. Note that the fused integral ()2.8|) admits an explicit 
evaluation as soon as the folded integral ()2.7|) admits some explicit evaluation. 
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Remark 2.1. A summation formula F = J2m&z n f( x + m ) f° r some function / : M n — »■ C, 
with F independent of x, implies the identity 

f(x) dx = F 

by the folding technique. This observation was extensively used by Gustafson [T^] to 
obtain evaluations of multivariate trigonometric integrals from known multivariate bilateral 
summation formulas. 

The simplest example of fusion and folding of integrals arises when (p(z) = 1 = ip{z). 
Then 

(z+m) 



Ei£±HLL q / s Si 
g 2 — § T [rz)q 2 



m=—oo 



is entire, hence 

(2.10) + (z) = CT_i(z) 

T 

is the corresponding sum ()2.6|) for some constant C. The integral (|2.7|) is 

C [ $_i{x)dx 



which equals the constant C in view of the explicit Fourier expansion ()2.5|) of the Jacobi 
theta function. The fused integral ()2.8|) thus becomes 



q^~dx, 

oo 



and the above analysis shows that it equals the constant C. The fused integral is up to a 
change of integration variable the Gauss integral, whose well known evaluation is given by 

q 2 dx 



yf—lT 

We thus conclude that the C = /_ iT m (|2.10|) . Returning to the underlying bilateral sum, 
we obtain 

oo 

\J—%T 
m=—oo v 

which is the well known inversion formula for the Jacobi theta function. The Jacobi 
inversion formula can be rewritten as 

2 

#_i(z) = V '—it q^ '§ t (tz), 

t 

or equivalently, with z replaced by z/t, 

i?_ 1(2/7-) = V—iTq 2 $ T (z). 
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Note that in this example, we used the explicit evaluation of the fused integral to derive 
the underlying bilateral summation formula (J2.11|) explicitly. Conversely, the evaluation of 
the Gauss integral follows from the Jacobi inversion formula and the elementary integral 



In subsequent sections we show how hyperbolic versions of beta integrals can be derived 
by the fusion procedure. In these cases, ip(x) = ip^ T ~i(x) is a slight modification of the 
(one-periodic variant) of the integrand of the trigonometric version of the beta integral with 
respect to base q, while (f)(x) is roughly ^(tx)" 1 (thus with respect to the base q); this is 
precisely what one would expect from the explicit formula (jl.7|) expressing the hyperbolic 
gamma function as a quotient of trigonometric gamma functions. 

3. The hyperbolic Ramanujan integral 

In this section we derive a hyperbolic analogue of the Ramanujan integral The 
derivations in this section are based on the (/-binomial theorem 



An elementary proof of the g-binomial theorem can e.g. be found in ^Hj. Of the many 
known proofs of the Jacobi triple product identity, we note that there are several which es- 
sentially only uses the g-binomium theorem, see e.g. £Q. We have chosen to give full details 
of all other identities we encounter since it clarifies the techniques leading to hyperbolic 
beta integrals. 

3.1. Ramanujan's summation formula. In this subsection we consider a nontriv- 
ial example of the bilateral sum construction ()2.6|) . It leads to an elementary proof of 
Ramanujan's {^i summation formula. This method of proving bilateral sum identities has 
been applied in several closely related setups, see e.g. 

Let t e EI and < \a\, \b\ < 1. We consider the --periodic entire function 



(2.12) 




[az; q) 

(3.1) i0 o (a;-;g,z) = / . \ > I?I,N<1, 




the Jacobi inversion formula (j2.11|) and the Jacobi triple product identity 
(3.2) tf T (z) = (q, -qa exp(2vuz), -q* exp(-2iTiz); q)^, \q\ < 1. 



(3.3) 4>{z; a, b) = <j>(z) = (-aq^ +z , -bq^ z ; q) 



and the associated entire function <\>[z\ a, b) = 4>(z) defined by 

~ £. 
<f>(z; a, b) = cj)(z; a, b) q 2 . 

The function satisfies (12.2)1 . as well as <\>{z + 1) = t(z)(/)(z) with 
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As explained in §2, we now force to become one periodic by considering the bilateral 
series 



(f) + (z)= Yl <f>( z + m ). 



which is easily seen to converge absolutely and uniformly on compacta of C to an entire 
function <p + (z) due to the conditions on a and b. Thus we conclude 

(f)+(z) = Cli(z) 

T 

for some constant C G C. To find an explicit expression for the constant C, we first express 
(f) + (z) in terms of the i^i bilateral series. Using <f>(z + 1) = t(z)<f>(z), we can rewrite 4> + {z) 
as 

oo 

(3.4) (f> + {z) = 4>(z + m) =t + ( y z)4>(z) 

m=— oo 

with 

oo m— 1 —1 —1 1 

<+ w = e n *(* + *) + e n jTTTi) 

l J - / m=0 fe=0 m=-ook=rn V ' 

= i^i -ag5+ 2 ; 6) 

(empty products equal 1 by convention). Now note that (for generic a and b), t + (z) can 
be evaluated at z G C satisfying q z = —qia^ 1 , in which case t + (z) reduces to a i0o series. 
By the g-binomial theorem we thus conclude 

t+(z)l («/«;«)«, 



On the other hand, 

0(2) 0(2) {-aq2 +z ,-bq2 z ;q) oQ 

where we used the Jacobi inversion formula and the Jacobi triple product identity for the 
last equality, hence 

"«•=-«*«-' (0.6; q ) m 

We thus conclude that 

(3.6) C= 1 (a6;9L 



V^F (0,6; g)„ 

The resulting bilateral summation formula 

(3.7) + (^;a,6)= ]T 0(z + m; a, 6) = V 
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is Ramanujan's 1^1 summation formula, written in the convenient form (|2.6J) . 

3.2. The Ramanujan integral. We consider in this subsection the natural trigonometric 
beta integral associated to Ramanujan's 1^1 sum 1)3. 7|) . 

Again we take r G H, so that |g| < 1. Let T = {z G C | |2| = 1} be the positively oriented 
unit circle in the complex plane. Choose parameters c, d G C satisfying 1 < \c\ < |g~ 5 | 
and < \d\ < \q~^\- We compute the trigonometric integral 

# T (z) ^ 1 f {q,-q^z~\-qh;q) oQ dz 



z 

' CO 



'o (— <pcexp(— 2nix), — q 2 dex.p(27rix); g) 27rz J t ^—q\ cz ~i ; —q\dz;q) 

by shrinking the radius of the integration contour T to zero while picking up residues 
(the second expression is obtained from the first by the Jacobi triple product identity and 
a change of integration variable). This computation is a simplified version of a general 
residual approach to basic contour integrals developed by Slater, see O §4.9] for details 
and references. 

By the conditions on the parameters, the poles — cq^ +m (m G Z>o) of the integrand lie 
inside T, while the poles — d~ 1 q~2~ m (m G Z> ) lie outside T. The condition |c| > 1 then 
ensures that 

1 [ (g, -q^z-\-qh; q) OQ dz _ Rgg / (g, -qh~\ -qh; g) c 

2ttz 7 t (-gl^- 1 , -g^dz; g)^ z m=Q 2= _ C( ^+™ \ (-gicz" 1 , -q^dz; q)^ z 
cf. [TBI §4.9]. The residues are easily computed, leading to 

1 f (q, ~q^z-\ -qh; q) OQ dz _ (gc, c~\ q) ^ ~ (gcci; g) m 



2vn 



(-glcz-^-gl^g)^ 2 (gcd; q) ^ ^ (g;g) 

(gc, gd; q) 



{qcd; g)^ 

by the g-binomial theorem. The conditions |c| > 1 and \d\ > may now be removed by 
analytic continuation. This proves the Ramanujan integral (|1.3jh 

(3.8) / _, -1W ix - ^' " d - ^ 



(— g2cexp(— 2-kix), — g2c?exp(27rix); g) (get/; g)^ 

for r G H and |c|, \d\ < \q~^\- The connection with the integrals from Ramanujan's lost 
notebook will become apparent in the following subsection. 

3.3. Fusion. The Ramanujan summation formula (J3.7j) and the trigonometric Ra- 
manujan integral ()3.8)) can be fused as follows. 

Proposition 3.1. Let r G H, < |a|, |6| < 1 and |c|, \d\ < \q~^\. Then 

11 x 2 

(3Q) f°° (-aq2 +x ,-bq2- x ;q) oD q S 2- ^ \ (ab-q)^ {qc.qd-q)^ 



(-g5cexp(-27r?x), -g^(iexp(27rix);g) oo V^rr (qcd; q) ^ (a,6;g) ( 
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Proof. The integrand J(x) of ()3.9|) can be written as 

4>(x; a, b) 



J(x) 



(— <pcexp(— 2mx\ —q^dexp(2TTix); 5) 
Note that the denominator is one-periodic. Hence folding the integral gives 

J(x) dx = f - — — — — dx 

Jo (— gacexp(— 2irix), —qidexp(2mx); q) 

V Jo ° 1 T dx 



\f—ir (a, b; q)^ J (— <pcexp(— 2nix), — q2dexp(2irix); q) 

where we used Ramanujan's summation formula (|3.7|) for the second equality and the 
trigonometric beta integral ()3.8j) for the last identity. □ 

The proposition is also valid for a = b = 0, in which case the resulting integral 
(3.10) /~ ^4 -^<b- 1 



(— g2cexp(— 2mx), — qzdexp(2mx); g) \J —ir ( < qcd;q') oo 

is the fusion of the Jacobi inversion formula (j2.11|) and the trigonometric Ramanujan 
integral ()3.8|) . The integral ()3.10|) is one of Ramanujan's trigonometric analogues of the 
beta integral from his lost notebook, see [3], ^3 Exerc. 6.15] and references therein. 
Ramanujan's second trigonometric analogue of the beta integral is the special case c = d = 
of the proposition, 



(-aq2 +x ,-bq*- x ;q) oo q^ dx 

-00 



(ab;q) ( 



>/=*F(a,6; q)^ 



This integral is the fusion of Ramanujan's 1^1 summation formula ()3.7j) and the elementary 
integral flZZTJ) . 

3.4. The hyperbolic Ramanujan integral. The integral ()3.9|) turns out to be more 
than just a particular (q, q )-bibasic extension of the Ramanujan integral. As we will see 
in this subsection, it is the crucial intermediate step towards the hyperbolic Ramanujan 
integral. 

We specialize the parameters in the fused Ramanujan integral (13.9)1 to 

(3.12) a = q a , b = q 13 , c = exp(— 2iria), d = exp(— 2m(3). 

We assume that r G H and that a, {3 G C are such that the corresponding parameters 
a, b, c, d satisfy the parameter constraints of Proposition ^. 11 The fused Ramanujan integral 
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()3.9|) can then be written entirely in terms of the r-shifted factorial (jl.7j) . 



(3 i3) / £ dx 1 ^i±3ii±M 

1 l ± + ± + a + x,± + ± + [3-x;T} 00 y/~FF [± + a + (3; r]* 



-oo L2 



By the reflection equation (J6.14)) we can rewrite f)3.13j) in a form which closely resembles 
the trigonometric Ramanujan integral (TH 



(3|J) l l| + h + gj| + E ~ ^ r](X) J X - q ' L + 7 + # T ] 



L2 + 27 + a + sc, g + 27 + P ~ x 'i r ]°o v 7 -^ 1 [7 + a + P\ T ]oo 

One may view (|3.14|) already as an hyperbolic analogue of the Ramanujan integral. It is 
though not completely satisfactory, because it does not extend to r G M<o (in which case 
the deformation parameters q and q satisfy |g| = \q \ = 1). To include this important case, 
we first need to rotate the integration cycle K. to — iR about the origin. 

For the rotation of the integration cycle it is convenient to restrict attention to r G 
EnC_, where C_ (respectively C+) is the open left (respectively right) half plane in C, 
and to parameters a, (3 G C satisfying 

(3.15) 0<a,(3<^. 

With these parameter constraints, the corresponding parameters a,b,c,d (see f)3.12jl ) au- 
tomatically satisfy the conditions of Proposition 13.11 We write 

\- + — + z - + — - z-t] 

j/ \ _ i2 T 2t T ^' 2 T 2t 1 l°° 



[\ + I + a + z,\ + I + (3 - z^U 

for the integrand of ()3.14|) . By the zero and pole locations of the r-shifted factorial (see 
(jSHIJ) and ()6.12|) respectively) we conclude that the poles of I(z) are contained in the two 
sets 

111^ 
-« + -- — + -Z< + Z> , 
(3.16) J 2r r 

P-~ + 7T + - Z >o + Z <o- 
Z Zr t 

By the conditions r G EI fl C_ and ()3.15|) . the first sequence —a + ~ — ^ + -Z< + Z> 
lies in the interior of the first quadrant of the complex plane (so in EI fl C+), while the 
second sequence (3 — \ + + ^Z> + Z< lies in the interior of the third quadrant of the 
complex plane (so in (— H) fl C_). Hence the integrand I(z) of 1)3.14)1 is analytic in the 
second and fourth quadrant of the complex plane. We now clockwise rotate the integration 
contour IR to — zR in (|3.14|) about the origin. We thus stay within the union of the second 
and fourth quadrant and do not pass poles of I(z). To rigorously verify that the integral 
evaluation does not alter, we need to take the asymptotic behaviour of the integrand I(z) 
into account. 
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Lemma 3.2. Suppose that r G Hn C_ and that the parameters satisfy (I3.15J1 . Then there 
exists a constant C G M>o such that 

\I(z)\ < C\qP z \ if Re{z) > and Im(z) < 0, 

\I(z)\ < C\q~ az \ if Re(z) < and Im(z) > 0. 

Proof. Using the reflection equation (j6.14j) we can write I(z) as 



I(z) = 



(3.17) 



[| + h + ot + z; 



C 2 q 



2t 
I 1 + ~ 



a 



z; r] c 



Z\ T 



for certain irrelevant ^-independent nonzero constants C\ and Ci- For the asymptotics in 
the fourth quadrant we use now the first expression for I(z) in ()3.17|) . Using the infinite 
product expression (jl.7j) for the r-shifted factorial, we may thus write 



(3.18) 



Hz) 



(-qhd 1 exp(-2^);g) oo (-ql +a q z ;q) oo pz 
— — ; — : — . ; ^ ; t~; ; — 1 



(-q 2 C exp(-27r^);g) oo (-ga V;?) c 

with c and d given by ()3.12j) . Now note that z i— > exp(— 27ri2;) maps the fourth quadrant 
G C | Re(z) > 0, Im(z) < 0} of the complex plane into the closed unit disc 



Since furthermore \q 2 
that 



D 

< 1 and |c| 



{w G 



\w\ < 1}. 

1 by the parameter constraints (|3.15jl . we conclude 
(-q ^d' 1 exp(-2iriz); q) 



{—q acexp(— 2niz); g) 

is uniformly bounded on the fourth quadrant of the complex plane. Observe that z 1— > q z 
maps the fourth quadrant of the complex plane into the closed unit disc D since r G HflC_. 
Furthermore, |<p -/3 | < 1 by the parameter constraints (J3.15)) . hence 

(-q^ +a q z ;q) c 



(-<P Pq z ;q) c 

is uniformly bounded on the fourth quadrant of the complex plane. In view of 1)3.18)1 this 
proves the uniform asymptotics for I (z) with z in the fourth quadrant of the complex 
plane. The asymptotics of I(z) in the second quadrant of the complex plane is determined 
similarly using the second expression for I(z) in (j3.17j) . □ 

We have the following direct consequence of Lemma 13.21 

Corollary 3.3. Let r G H fl C_ and suppose that the parameters a,/3 G C satisfy (|3. 15)1 . 
Let < e < 1 be the growth coefficient 

e = exp(— 27r min(a, (3) m) 
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with m the strictly positive constant 

m = min (lm(Te t6 )) . 

0e[-f,o] 

Then there exists a constant C G R>o such that 

\I(z)\ < Ce^ 

for z G C satisfying Re(z) > and Im{z) < 0, as we// as /or z G C satisfying Re(z) < 
ana 1 Im(z) > 0. 

We keep the assumptions r G HflC_ and (|3.15|1 . In view of Corollary 13. 31 we may apply 
Caucliy's theorem to rotate clockwise the integration contour R in f!3.14|) to — iR about 
the origin without altering its evaluation. Relaxing the parameter constraints leads now 
to the following main result. 

Theorem 3.4 (Hyperbolic Ramanuj an integral). Letr G C with Re(r) < and Im{r) > 0, 
and suppose that the parameters a, (3 G C satisfy 

(3.19) Re(ra),Re(T/3) < 0, Re(a - - + — ), Re{(3 - \ + -*-) < 0. 
Then 

(3.20) r 1 \ + ^ + ^\ + h- z ^]oo dz g-ig"[J + a,J + /3;r] 00 



L2 + 27 + a + ^, 5 + 27 + /3-£;r]oc " [;; + « + t] 



Proof. Let r G HI fl C_. We already argued that (J3.20p is valid for parameters a and 
/3 satisfying the constraints (|3.15j) . By analytic continuation one easily verifies that the 
integral evaluation (J3.2Uj) is valid under the milder constraints (J3.19j) . In fact, in view 
of ()6.15|) and ()6.16|) the requirements that the integrand I(z) decays exponentially for 
z — > ±ioo and that the sequence of poles —a + | — j- + ^Z< + Z> (respectively (3 — 
| + 27 + 7^>o + Z< ) of I(z) is contained in C + (respectively C_), lead to the conditions 
f!3.19|) on the parameters a and f3. 

Let t' G R<o and choose a, /3 G C satisfying ()3.19|) with r = r'. Let < 5 < | and 
define the compact subset 

K= {r(0)|0G [0,5]} C C_ 

with t(4>) = t' exp(— icf)). By choosing 5 small enough we may assume that the parameter 
constraints (j3.19j) are valid for all r G K. Observe that K n HI = {r(0) | < < 5} 
and that K fl R<o = r(0) = r'. By the reflection equation (|6.14|) and by the asymptotic 
estimates ()6.15|) and (j6.16|) for the r-shifted factorial, the weight function I(z) in (|3.20|) 
has the asymptotic behaviour 

(9(exp(— 2-Kirctz)) , if lm(z) — > oo, 

(9(exp(27rir/?2;)) , if Im(^) — > — oo, 

uniformly for t E K and for Re(z) in compacta of R. Consequently, 

|J(zx)| = 0(e N ), x^±oo 
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uniformly for r G K, with growth exponent 

e = exp 27r max (Re(ra),Re(r/?)) ) < 1 



rel< 

depending only on a, (3 and K. Thus if we take r = r(0) = KnM. in (|3.20p . then Lebesgue's 
dominated convergence theorem allows us to interchange the limit <fi \ with the integral 
over iR. This yields (ET27IJ) for r = r'eR <0 . □ 

4. The hyperbolic Askey- Wilson integral 

In this section we apply the folding and fusion techniques to derive an hyperbolic ana- 
logue of the Askey- Wilson integral. For proofs of the underlying summation formula and 
the underlying trigonometric Askey- Wilson integral (|1-4|) we refer now directly to the lit- 
erature (in practice we often refer to the book ^Hj of Gasper and Rahman, in which one 
can find further detailed references to the literature). 

4.1. Bailey's summation formula. We take r 6 H, so that |g| < 1. For the Askey- 
Wilson integral and its hyperbolic analogue, the role of the function in the underlying 
summation formula (see §2) is 

( 41 ) UU( s ^ ±Z ^)oo 



(q^^-q^.q^-q)^ 
As in §2, we define 

~ id 
(j)(z) = <f>(z) q 2 . 

We assume for the remainder of this subsection that the parameters Sj are nonzero complex 
parameters and satisfy \q~ 3 SiS2S3S4\ < 1. A direct computation shows that 

$(z+l)=t(z)j(z) 

with 

( ) = giWj (1 ~ 1 l+z ) (1 + Q 1+z ) A (1 - * 3 - V + 

<f (i-g*)(i + g*) f}, (i-a jq *) 

Due to the conditions on the parameters the bilateral sum 

oo m— 1 —1 —1 

m=0 k=0 m=— oo k=m 

fq 1+z ,-q 1+z ,s^q 1+z ,s 2 l q 1+z ,Ss l q 1+z ,sfq 1+z s 1 s 2 s 3 s i \ 



! '* ' q z ,-q z ,s 1 q z ,s 2 q z ,s 3 q z J s 4 q z ' g ' g3 y 

converges absolutely and uniformly on compacta away from the Z< -translates of the poles 
and the Z >0 -translates of the zeros of t(z), and it extends to a meromorphic function in 
z6C. Hence the bilateral sum 



4> + {z) = (f){z + m) 
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converges to a meromorphic function in z G C, given explicitly by 

<f> + (z) = t + (z)j(z). 

We now take a shortcut compared to our analysis of the Ramanujan integrals by directly 
applying Bailey's summation formula jTHJ (5.3.1)] for the explicit evaluation of t + (z), 

t+( z ) = (?L q )°° ni<*<m<4(9 _1 ^sm; q)^ (g 1±2z ; q)^ 



(q 3 s 1 s 2 s 3 s^q) oo UU( s ^ ±z '^)oo' 

Since 

(4.2) (u 2 ; q) oo = (u, -it, <^u, g)^, 

the meromorphic function <f) + (z) = t + (z)4>(z) can now be written as 

(4.3) <p + ( Z ) = ^n^^wi ^ i{z) 

y/-ir {q~ 3 s 1 s 2 s 3 s^q) oo 

after a straightforward computation using the Jacobi triple product identity and the Jacobi 
inversion formula. Written in the form ()2.6j) . we thus have <p + (z) = &(z)'&_i(z) with the 

T 

elliptic function <&(z) being the constant 

1 IlKk^i?" 1 ^^;?)™ 



(4.4) $ 



\T i vt (q 3 s 1 s 2 s 3 s 4 ;q) ( 



4.2. Fusion. In this subsection we fuse the trigonometric Askey- Wilson integral (jl.4j) with 
Bailey's summation formula ()4.3j) . For proofs of the trigonometric Askey- Wilson inte- 
gral (jl.4|) . see [7], [T3J §6] and references therein. 

Proposition 4.1. Let r G HI and Sj,t k G C (J, k = 1, . . . , 4) wift |g _3 SiS2S3S4| < 1 and 
\t k \ < 1 (jfe= 1,...,4). Then 

(exp(±2irix)i-exp(±2-Kix),g^ exp(±2mx); q ) oo ^ (sjg ±3! ; g)^ 

2 {tit 2 t 3 t 4 ;q) ^ j_j (g -1 ^™;?)^ 



v^(g 3 sis 2 s 3 s4;g) 00 1 < fc A < i< 4 

Proof. The integrand «/(x) can be written as 

J(x) = 0(x)w(x) 
with </>(x) as in §4.1 and with w(x) the one-periodic function 

(exp(±27rzx), — exp(±27rzx), q 2 exp(±27ux); g) 



n i= i(^exp(±27rix);g) oo 

Note that J(x) is regular on R since the real zeros of the factor (q 1±x , q^ ±x ; g) in the de- 
nominator of J(x) are compensated by zeros of the factor (exp(±27rza;), — exp(±27ro); q ) 
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in the numerator of J(x). Thus folding the integral and using Bailey's sum (14. 3j) . we 
obtain with $ the constant (J4.4j) . 



J(x) dx = J (f)(x)w(x) dx 

4> + (x)w(x) dx = $ / $_i(x)w(x) dx 



by Fubini's theorem. Using now the Jacobi triple product identity and (|4.2jl . the resulting 
integral over the period cycle [0, 1] is essentially the trigonometric Askey- Wilson integral 
(II- 4(1 in base q. Its evaluation gives the desired result. □ 

Remark 4.2. Under the parameter constraints of Proposition 14. II we have 

f(x)J(x)dx = & / f (x)d _i(x)w{x) dx 



for a regular one-period function / : R — > C, where we used the notations introduced in 
the proof of Proposition 14.11 The latter integral is essentially the integral of / against the 
Askey- Wilson measure depending on the four Askey- Wilson parameters tj (j = 1, ... ,4) 
and with base q. In particular, the Askey- Wilson polynomials [7j, [T3J §7.5] 



q n ,q n H^tst^, ti exp(27iix), ti exp(— 2nix) _ 

t\t2, tits, titi 



Pn\ x ) — 403 ++++++ ) 9) 9 ) n ^ ^>0 



form a basis of the space C[cos(27rx)] of polynomials in cos(27rx), and they are orthogonal 
with respect to the complex pairing 

/oo 
Pi(x)p 2 (x)J(x)dx. 
-oo 

The "quadratic norms" (p n ,Pn) {n £ ^>o) can be explicitly evaluated using the quadratic 
norm evaluations of the Askey- Wilson polynomials, see [Zj and [TBI §7.5] (the special case 
n = is the statement of Proposition I4.1J1 . 

The proposition is also valid for (s%, s 2 , s 3 , s 4 ) = (q, —q, q^ , 0), in which case the resulting 
integral 

(exp(±27ra;), — exp(±27rix), q 2 exp(±27ra;); q)^ 
U 4 j =i(t j exp(±27rix);q) oo 

2 {tit 2 t z t A ;q) 



v /z ^ I ni<ife<m<4(^m;g r ) E 



is the fusion of the Jacobi inversion formula (j2.11j) and the trigonometric Askey- Wilson 
integral JO}. This integral may be viewed as a one- variable Macdonald-Mehta integral, 
see and jS]. In particular, the special case (ti, t 2 , t 3 , t 4 ) = (0, g 2 y q k ; — g fc ) with fc G 



>Q 
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leads to Cherednik's [H] one-variable Macdonald-Mehta integral 

(exp(±4mx);q 2 ) j 2 {q 2k ;q 2 ) ( 

q 2 ax 



, (q 2k exp(±ATTix);q 2 ) oQ ^f^rr (q Ak ; q 2 ) ^ 

In a similar way we would like to put (t\, t 2 , £3, t 4 ) = (1, — l,q?,0) in the fused Askey- 
Wilson integral (Proposition 14. If) . but this is not allowed since the integrand is irregular 
on the real line for these parameter values. A suitable regularization of this integral yields 
an integral evaluation which is closely related to a trigonometric beta integral of Askey 
5, §3], and which has a similar appearance as Etingof's §3.5] Macdonald-Mehta type 
integral with real integration cycle. The result is as follows. 

Corollary 4.3. Let r G HI and Sj G C (j = 1, . . . , 4) with \q~ 3 SiS2S3Si\ < 1. Then 

lim r +0 ° UUi^^oo fi dz i m^n^Or 1 ^;^ 



c\fl Jie-oo (g 1±2 ,-g 1±z ,g5 ±z ;g) oo 4^vr (g 3 sis 2 s 3 s 4 ;g) oo 

Proof. The integrand <j){z) is regular on ie + R for e G M>o sufficiently small. Folding the 
integral yields 

»ie+oo />1 />1 

[z) dz = / </> + (ze + x)c?a; = $ / $_i(ie + x)dx 



o 



by JO} and (jPjh By (gH) we have 

-1 /•! 



lim / $ i (ie + x)c?x = / $ i(x)<ix=l, 
'Wo " Jo " V 

which completes the proof. □ 

4.3. The hyperbolic Askey- Wilson integral. In this section we derive the hyperbolic 
analogue of the Askey- Wilson integral and, as a special case, the analogue of the Askey- 
Wilson integral for |g| = 1. The techniques are the same as for the derivation of the 
hyperbolic Ramanujan integral in §3.4. 

For the moment we assume that r 6 Hfl C_. As in the previous section we write 
q = q T = exp(2mr) and q = q T = exp(— 2m/ t), and q u = exp(2mru), q u = exp(—2mu/r) 
for u G C. We fix four parameters Tj (j = 1, . . . , 4) satisfying the parameter constraints 

l-r,-GHnC + , (l-Tjjrei (j = l,...,4), 
(4-5) , 

(ri + r 2 + T3 + T4 — 3 — r )r eHnC_. 

If we take parameters Tj G C satisfying 

Re(r i ) < 1, Im(rj) < 0, Re(n + r 2 + r 3 + r 4 ) > 3, 

then sufficient additional conditions on their imaginary parts to satisfy the parameter 
constraints (|4.5j) are 

(4.6) Imfo) > (l-Re( Ti ))^4, Im(n +r 2 +r 3 +r 4 ) > R e (r 1 +r 2 +r3+r 4 -3) R< ( T) 



Re(r) ' v 1 z a * y v 1 z a * ; Im(r) 
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The existence of parameters satisfying these inequalities follows from the fact that the right 
hand sides of the inequalities in (J4.6|) are strictly negative. 

For fixed parameters Tj satisfying the parameter constraints (j4.5j) we associate eight 
parameters Sj, tj (j — 1, . . . , 4) by 

(4.7) Sj = q T \ tj = exp(-27rzr i ), (j = 1, . . . , 4). 

These eight parameters satisfy the parameter requirements of Proposition 14.11 The result- 
ing fused Askey- Wilson integral can then be written in terms of the r-shifted factorial (jl.7)l 

as 

(a a\ f°° [1±x,\±x,1 + ^±x;t] oc ^ 2 [n + r 2 + r 3 + r 4 - 3; r]^ 

v 4 -"/ / 7^1 — r ; 1 q 2 ax 



Y[j =1 [Tj ± x; r]oo V^rUi<k<m<A r k + r m - 1; r] 

Using the reflection equation (j6.14|) and 

/ x rl 11 1 

(4.9) - ± x, 1 ± x, - + — ± x, 1 + — ± x; t}^ = [1 ± 2x; t m 

2 2 2t 2t 

(which follows easily from (jl.7|) and (J4.2j) ). we can rewrite (|4.8|) as 
(41Q) [l±2x;r] 00 dx _2q~™q™ [n + r 2 + r 3 + r 4 - 3; r]^ 



U j=1 [Tj ±X;r]oo ^f^VT Ul<k<m<ii T k + T m-l;r}^ 

which now has a very similar appearance as the trigonometric Askey- Wilson integral (jl.4|) . 
The numerator of the integrand 

(4.1D nz) = ^ff 1 -, 

can be simplified as follows. 
Lemma 4.4. 

[1 ± 2z; r]oo = -iq^q~^q 2z2 (q~ z - g 2 ) (exp(27rzz) - exp(-27T22:)) . 
Proof. By the reflection equation (|6.14j) we can write 



i~i 1 ( 1 1 2A 2 \2z + 1:t 
[l±2z-T} OQ = q-^q^q^-~ 2z > 1 \ 

[2z + r _i ; t 



OO 

) ' ]oo 



Now apply the functional equations ()6.10|) to the numerator and denominator to get rid of 
the r-shifted factorials. Simplifying the resulting elementary expression gives the desired 
result. □ 

By the zero and pole locations of the r-shifted factorial (see (|6.11|) and (j6.12J) respec- 
tively) we conclude that the poles of I(z) are contained in the eight discrete sets 

(4.12) ± Uj + iz> + Z <0 ) (j = 1, ... ,4). 

By the parameter constraints r G HI fl C_ and ()4.5|) . the four sets in ()4.12|) with plus sign 
are contained in the interior of the third quadrant of the complex plane, and consequently 
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the four sets with minus sign are contained in the interior of the first quadrant of the 
complex plane. 

Lemma 4.5. Suppose that r G In C_ and that the parameters Tj (j = 1, ... ,4) satisfy 
(14.5)1 . Then there exists a constant C G R>o such that 

y ( 1~1 + T~2 + T3 + T4 — 3 — T ~ 1 ) Z 



\I(z)\ < C 



Q 



for z G C satisfying Re(z) > and Im(z) < 0. 

Proof. By Lemma f4. 41 the reflection equation ()6.14|) and by (|1.7|) we can rewrite the inte- 
grand I(z) as 

jf z \ =CT\ (^ exp ^ 27rir ^ exp ^ 27r ^ ; ^)oo (g^gj_g)oo 
7=i ( ex P(- 27rir i) ex P(- 27r ^);5') 00 (q l ~ Tj q z ; q)^ 

x (1 - exp(-4vr^)) (l - g 2z ) g (n+^+7 3 +r 4 -3-r- 1 ) 2 

for some irrelevant nonzero constant C. The proof is now analogous to the proof of Lemma 
13.2) since | exp(— 2mTj) \ < 1 and |g 1 ~ T -'| < 1 by the conditions (J4.5j) on the parameters. □ 

Note that a bound for the weight function I(z) for z G C satisfying Re(z) < and 
Im(z) > can be immediately deduced from Lemma l4"31 using I(—z) = I(z). The following 
corollary follows from Lemma 14.51 and the conditions ()4.5|) on the parameters. 

Corollary 4.6. Let rein C_ and suppose that the parameters Tj satisfy the parameter 
constraints (|4.5|) . Let < e < 1 be the growth coefficient e = exp(— 2nm) with m the 
strictly positive constant 

m = min (lm((ri + r 2 + Tq + r 4 — 3 — r^Mre 16 )) . 
0e[-f ,o] v v /y 

Then there exists a constant C G M >0 suc/j t/iat 

\I{z)\ < Ce lzl 

for z G C satisfying Re(z) > and imfz) < 0, as well as for z G C satisfying Re(z) < 
and Im(z) > 0. 

We keep the assumptions r G EI fl C_ and f)4.5|) . In view of Corollary 14.61 we may apply 
Cauchy's theorem to rotate clockwise the integration contour R in (14.10)1 to — iM. about 
the origin without altering its evaluation. Relaxing the parameter constraints leads now 
to the following result. 

Theorem 4.7 (Hyperbolic Askey- Wilson integral). Let r G C with Re(r) < and 7m(r) > 
0, and suppose that the parameters Tj G C (j = 1, . . . , 4) satisfy 

(4.13) Re(Tj)<l (j = l,...,4), J Re((r 1 + r 2 + r 3 + r4-3)r) < 1. 

Then 

(4M) f l °° [l±2zyrU dz _ 2g~AgA [ Tl + r 2 + r 3 + r 4 - 3; r]^ 

J-ioo Y\ A j=1 [Tj ± 25 T~]oo V 7 -^" ril<fc< m <4[ r fe + r m - I! t]oo 
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Proof. For r G in C_ we already argued that the integral evaluation f)4. 14j) is valid for 
parameters Tj satisfying the parameter constraints (J4.5|) . By analytic continuation one 
proves that the integral evaluation 1)4.14)1 is valid for parameters Tj satisfying the milder 
constraints (J4.1H)) . In fact, by (j6.15j) and (jfi.lfij) the open convex parameter domain of 
parameters (t\, T2, T3, T4) satisfying the constraints ()4. 13|) is the maximal extension of the 
parameter domain ()4.5|) such that the associated integrand I(z) decays exponentially as 
z — > ±ioc and such that the sequence of poles Tj + ^Z> + Z <0 (j = 1, . . . , 4) lie in C_. 

The extension of the integral evaluation (|4.14j) to r G M<o follows now in a similar manner 
as for the hyperbolic Ramanujan integral (see the proof of Theorem l3.4|) using the reflection 
equation (j6.14|) and the asymptotics (j6.15J) and (j6.16J) of the r-shifted factorial. □ 

Remark 4.8. The Askey- Wilson polynomials (see Remark |4.2|) do not remain an orthogonal 
system in the hyperbolic case. Indeed, the Askey- Wilson polynomials are one-periodic 
hence bounded on R, but they grow exponentially in the imaginary direction (with growth 
rate exp(27rn|Im(z)|) where n is the degree of the Askey- Wilson polynomial). The rotation 
of the integration cycle K to — z'R in the complex orthogonality relations of Remark 14.21 
is only allowed if the integrand has uniform exponential decay in the second and fourth 
quadrant of the complex plane. For given parameters satisfying (|4.5|) . this will thus only 
be the case for certain low degree Askey- Wilson polynomials. 

On the other hand, in |2E] an hyperbolic analogue of the Askey- Wilson function is 
defined, which is expected to play the role of integral kernel for some generalized q- Fourier 
transform with \q\ = 1. 



5.1. A weak g^g summation formula. We take r G H, so that \q\ < 1, and we fix 

five generic nonzero complex parameters tj (j = 0, ...,4). It is convenient to write 
A = to^i^2^3^4 for the product of the five parameters. For the hyperbolic analogue of 
the Nassrallah-Rahman integral, the role of the function (j) in the underlying summation 
formula (see §2) is 



5. The hyperbolic Nassrallah-Rahman integral 



(5.1) 



4>(z) 



q 1±e , q^ ±z , Aq~ 4±z ; q) 



00 



As in §2, we define 



z) = (j)(z) q 2 . 



A direct computation shows that 



4>(z + l)=t(z)(f>(z) 



with t(z) given by 



t(z) = q 



(1 



- q 1+z ) (1 + q 1+z ) (1 - Ar 4 + 2 ) A (1 ~ ^V +Z ) 
(l ~ Q z ) (l + Q z ) (l _ A~ 1 q 5+Z ) 11 (l-tj<f) 
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The bilateral sum 

oo m—1 —1 —1 

m=0 k=0 m=— oo k=m 

fq 1+z , -q 1+z , t Q l q 1+z , t^q 1+z , t^ l q 1+z , t^q 1+z , t^q 1+z , Aq~^\ 
8 8 V q z ,-q z Mq z Mq z Mq z Mq z Mq z ,A- l q^ z ' y ' q 

converges absolutely and uniformly on compacta away from the Z<o-translates of the poles 
and the Z>o-translates of the zeros of t(z), and it extends to a meromorphic function in 
z6C. Hence the bilateral sum 



4> + {z) = cj)(z + m) 



converges to a meromorphic function in z e C, given explicitly by 

<P + (z)=t + (z)4>(z). 

Unfortunately, there is no explicit summation formula on the level. Instead we apply 
a formula that expresses t + (z) as a combination of two one-sided sums which only in- 
volve ^-independent summands. This can be done by applying [T3*l (5.6.2)] with its seven 
parameters a,b,c, . . . , g specialized to 

(a, b, c, d, e, /, g) = (q z , qt^ 1 , gt 2 \ qt^ 1 , qt^ 1 , qt \ q~ 4 A) . 

After a straightforward but tedious computation using ([4.2)1 . the Jacobi triple product 
identity and the Jacobi inversion formula, we then arrive at the explicit expression 

f (t q ±z , tn Y ±z ; q) } ~ 

with the ^-independent constants C\ and Gi given by 



1 lli=i(9 1 totj,qt 1 t j ;q) ( 



(5.2) {q- 5 At , q-*At- Q \ q% 2 ; q)^ 

x s W 7 (q 2 t 2 ; g 2 i V, g 2 to \ g 2 t V, ^o^ 1 . 9~ 3 ^o *l ?) 



and 



1 Ui^A-Hj^^ql 

Lo 



(5.3) (g 5 ^- V. (r 8 ^ 1 , t 7 ^ 42 ; ?L 

x 8 W 7 (q- 8 A 2 ; q^At^ 1 , q~ 3 At^, q- 3 At 2 \ q^At^ 1 , q^At^ 1 ; q, q). 
The meromorphic function <f) + (z) can thus be written as 

(t q ±z , to x q 1±z ; q) 

(5.4) 0+ * = $(z)'&_i(z), Hz) =C 1 + ( \°\ ± l ° \ ^ C 2 . 

{Aq-^^A^q^-q)^ 
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The function <&(z) is an elliptic function in z with respect to the periods 1 and — -. For 
later purposes it is convenient to rewrite 3>(z) as follows. Fix Tj G C (j = 0, . . . , 4) such 
that 

tj = q T > (j = 0,...,4) 
and denote a = r + T\ + r 2 + r 3 + r 4 , so that A = q a . We define new nonzero complex 
parameters tj (j = 0, . . . , 4) by 

tj = exp(-27rzr i ) (j = 0, . . . , 4), 

and we write A = tot^t^t^ = exp(— 2nia). Then 

(to exr>(±2niz), qtT, 1 exp(±27rzz): g) 

$(z) = d+ L PV C 2 g 7 , 

^ ^ (A exp(±2niz),qA- 1 exp(±2iTiz);q) oo 

9 1 \ 2 / 1 1 ' 



by a direct computation using the Jacobi triple product identity and the Jacobi inversion 
formula. 

5.2. Fusion. In this subsection we fuse the weak 8 \l/ 8 summation formula (|5.4|) with the 
Nassrallah- Rahman integral (|1.5p . For proofs of the Nassrallah- Rahman integral ()1.5|) 
itself, see [TH], [H3 (6-4.1)] and references therein. 

We fix r G i and generic complex parameters Tj. As in the previous subsection, we 
write 

(5.6) tj = q Tj , % = exp(-27rir i ) (j = 0,...,4) 

and A = q a , A = exp(— 27ria) with a = r + 7~i + r 2 + r 3 + r 4 . We assume throughout this 
subsection that Im(r,) < for j = 0, . . . , 4 and that Im(a) > Im(r _1 ). These conditions 
imply that 

|5|<1 (j=0,...,4), \A\>\q\. 
By ([4.9)1 . the reflection equation (j6.14j) and the definition (jl.7|) of the r-shifted factorial, 
we can write 



[l±2z,a-4±z;r} 



i i 



(5.7) — — ■ j-!— = q ^q^(j)(z)v(z) 

with <p(z) as in §5.1, and with v(z) the one-periodic function 

( 5 g) ^ _ ( ex P( ±2?r ^)> ~ exp(±27riz), q ? exp(±27rcz), A exp(±2niz);q) ^ 

n j 4 =oS ex p( ±27r ^);5')oo 

Note that (|5.7|) is regular on R in view of the conditions on the parameters Tj. By Fubini's 
theorem, the integral 

Z" 00 fl ± 2x, a - 4 ± x; rloo , Z 100 ~ , . . , 

(5.9) / 1 — ax = q ^q^ / <p{x)v[x) ax 

J -co Y\.j=wTj i -^i T ]oo </-oo 
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can now be folded, 

(j)(x)v(x) dx — / 4> + (x)v(x) dx 
Jo 

= C\ I v{x)d_\{x)dx 
Jo 

+ cw /'«(«) (l° ex p( ±2A )-5g' ex p( ±2 '")-?L tf _ l(j)ttt 

Jo (Aexp(±27iix),q A' 1 exp(±2nix);q) 

by ()5.4|) and (|5.5|) . Applying the Jacobi triple product identity and (|4.2jl . the remaining two 
integrals over the period interval [0, 1] can be evaluated by the trigonometric Nassrallah- 
Rahman integral (jl.5j) . This gives 

(5.10) / (j){x)v{x) dx= 1 V J J 1°° <^ d + C 2 g 7 TT L J Y~ 3 J ^ g 

n <fe<m<4(^m;5') oo { f = \(At J \qA-H j -q) c 

with the constants C\ and C*2 given by (|5.2|l and (|5.3|l . respectively. The following key 
lemma gives the evaluation of the remaining sum. 

Lemma 5.1. 



Ci + (W ] I 



Proof. Using the Jacobi triple product identity and the Jacobi inversion formula we can 
write 

\} 1 {Atj\qA-H j - 1 q) oo U^AtjWA-H^ 
with <5 given explicitly by 

Writing out the squares for 7 ()5.5|) and 5 gives 



7 = At- 1 + 20 -87-0 + ^(7-2- (9 + t- 1 )tj) +2 ^ r fe r m = -5, 

j=l 0<fc<m<4 



hence 



(5.11) S^SiL =Ci+C2 tt (£Mij^Y4 



26 JASPER V. STOKMAN 

with C\ and C 2 given by (|5.2|) and (|5.3jl . respectively. The right hand side of f!5.11|) can be 
evaluated by Bailey's summation formula (2-11-7)] with the six parameters a,b, . . . , f 
specialized to 

(a, b, c, d, e, /) = {q% 2 , q~*MZ \ q%\\ q%H 2 \ q%H^\ qW). 

This yields the desired result. □ 

Combining the lemma with ()5.7|) and ()5.10j) immediately implies the following result. 

Proposition 5.2. Let r G H and tj G C with Im(Tj) < (j = 0, ...,4) and Jm(a) > 
Im{r~ l ), where a = r + ti + r 2 + r 3 + r 4 . T/iere 

(5 12) Z" 00 [l±2s,a-4±s; rjoo ^ _ 2g"^g H YlUoi a ~ T i ~ ^ T U 



Yl j=0 [Tj ±x;t] oq y^vr Y[ G < k< m<A T k + T m ~ ^] t}^ 

Remark 5.3. The weak 8 \I/ 8 summation formula (|5.4jl involves an elliptic function $(z) 
depending nontrivially on the parameters tj. Consequently, the fusion of the weak 8 \l/ 8 
summation formula with the trigonometric Nassrallah-Rahman integral (|1.5|) only has an 
explicit evaluation when the parameters in the trigonometric Nassrallah-Rahman integral 
are matched to the parameters in the weak 8 \l/ 8 summation formula via (|5.6|) . This is an 
essential difference with the fused Ramanujan integral in §3.3 and the fused Askey- Wilson 
integral in §4.2. 

It is possible to fuse the trigonometric Nassrallah-Rahman integral (jl.5j) in base q with 
the Jacobi inversion formula (|2.11jl . yielding the one- variable Macdonald-Mehta type inte- 
gral 

(exp(±27rzx), — exp(±27rzx), q I exp(±27rza;), Aexp(±2nix); <f) ( 



(5.13) 



n i =o( t i ex p( ±27ria; )^)oo 

2 m=oK^?) ( 



q 2 dx 



UoKkKm^i^m; q) ^ 

with A = tot^tsti and with parameters tj G C satisfying \tj\ < 1 (j = 0, . . . , 4). 

5.3. The hyperbolic Nassrallah-Rahman integral. Let rein C_ with Re(r _1 ) < 
—1 and fix five parameters Tj (j = 0, . . . , 4) satisfying 

(5.14) 1 — Tj, a — 4 — t -1 Glfl C+, (1 — Tj)r, (a — 4)r G H 

for j = 0, . . . , 4, where a = r + Ti + r 2 + r 3 + r 4 . If the real parts of the parameters Tj 
satisfy 

(5.15) Re(r j ) < 1 (j = 0, . . . , 4), Re(a) > 4, 
then they satisfy ()5.14|) when the imaginary parts of Tj satisfy 

(5.16) (1 - Re(Tj)) ^ < Im(r,) < 0, k(fl)^ < -1 - Re^ 1 ) 
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for j = 0, . . . ,4. Note that the condition Re(r _1 ) < —1 implies —1 — Re(r _1 ) > 0, hence 
there exist parameters Tj satisfying (|5.15|) and (J5.16)) . 

The constraints 1)5.14)1 imply the parameter conditions of Proposition E21 We now rotate 
the integration cycle K. clockwise to — z'K. about the origin in the fused Nassrallah-Rahman 
integral ()5.12)) . By the zero and pole locations of the r-shifted factorial (see (j6.11)) and 
(16.12)1 respectively) and by Lemma f4. 41 the poles of the integrand 

(5.17) I{z)= [l±2z,a-A± Z] rU 



of the fused Nassrallah-Rahman integral (J5.12)) are contained in the union of the twelve 
discrete sets 

(5.18) ± (rj + ^Z> + Z <0 ^ (j = 0, . . . , 4), ± (-a + ^Z >0 + Z< 4 

By the parameter constraints r G HI fl C_, Re(r _1 ) < — 1 and (J5.14)) . the six sets in ()5.18j) 
with plus sign are contained in the interior of the third quadrant of the complex plane, and 
consequently the six sets in ()5.18j) with minus sign are contained in the interior of the first 
quadrant of the complex plane. Furthermore, the integrand I(z) is even, J(— z) = I(z), 
and by ()6.14|) and p. 7)1 it can be rewritten as 

on/ \ (exp(— 2nia) exp(— 2niz); q) (q 5 ~ a q z :q) 

I(z) = C (1 - q 2z ) (1 - exp(-47r^)) V V x ' x 



x 



(q exp(27rm)exp(-2vr^);g) oo (g a " 4 g 2 ;g) oo 

* (gex P (27rir 3 -)exp(-27r»z);g) oo (q* q z ; q) „ \ qZ exp{ _ 2mz) 
f} (exp(-27rzr i ) exp(-2niz); q) ^ (q^q*; q)^ 



for some nonzero constant C . With similar arguments as in the proof of Lemma 13.21 and 
Corollary 13.31 for Lemma f4. 51 and Corollary 14.6)) we conclude that for some C\ G M>o, 

|/(^)| < c ie i-i 

for z G C satisfying Re(z) > and Im(z) < as well as for z G C satisfying Re(z) < and 
Im(z) > 0, with growth exponent < e = exp(— 2i\vii) < 1, 

m = min ( (r - l)e %0 ) > 0. 

Thus we may apply Cauchy's theorem to rotate clockwise the integration cycle IR in ()5.12j) 
to — zR about the origin without altering its evaluation. Relaxing the parameter constraints 
leads to the following result. 

Theorem 5.4 (Hyperbolic Nassrallah-Rahman integral). Let r G C with Re(r) < and 
7m(r) > 0, and suppose that the parameters Tj G C (J = 0, . . . , 4) satisfy 

(5.19) Re(Tj) < 1 (j = 0,...,4), Re(a - r" 1 ) > 4, 
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where a = r + T\ + r 2 + r 3 + r 4 . Then 

(5>20) p [l±2z,a-4±z;T] 00 ffc _ 2q-*q* Il}^[a- ^ - 3;r] 



rij=o[ r i ± ^ T ]oo V 7 -^ Ilo<fc<m<4[ r fc + T m - 1; 1~]oo 

Proof. For r G HflC_ satisfying Re(r~ 1 ) < —1 we have already argued that ()5.20|) is valid 
for parameters Tj satisfying the constraints (|5.14|) . By analytic continuation we conclude 
that ()5.20|) is valid under the milder parameter constraints ([5. 19)1 . using similar arguments 
as in the proof of Theorem 13.41 

The condition Re(r _1 ) < — 1 can be removed by analytic continuation using the asym- 
potics (uniform for r in compacta) of the r-shifted factorial, see (jbM5j) and (jbM6j) . The 
extension of (|5.2()|) tor G M<o then follows as in the proof of Theorem 13.41 □ 

Remark 5.5. The hyperbolic Askey- Wilson integral (J4.14j) is formally a limit case of the 
hyperbolic Nassrallah-Rahman integral. For example, let r G M <0 and choose parameters 
Tj G C (j = 1, . . . , 4) satisfying (j4.13|) . Concretely, the r/s thus satisfy 

Re(rj) < 1, Re(ri + r 2 + r 3 + r 4 ) > 3 + r' 1 . 

Choose Tq G 1R<i such that 

t'u + Re(n + r 2 + r 3 + r 4 ) > 4 + r -1 , 

then any five tuple (r , ri, r 2 , r 3 , r 4 ) with Re(r ) = Tq satisfies the parameter constraints 
(I5.19J1 . We can now formally take the limit Im(ro) — > — oo in the hyperbolic Nassrallah- 
Rahman integral ()5.20|) while keeping Re(ro) = Tq fixed. By (|6.17|) . this formal limit of the 
hyperbolic Nassrallah-Rahman integral is the hyperbolic Askey- Wilson integral (|4.14|) . 

5.4. The hyperbolic degeneration of the elliptic Nassrallah-Rahman integral. 

We introduce the short hand notation t = (to,ti,t 2 ,t 3 ,t^). With the notations of §1, we 
write 

5 -21 A(z;t; Pl ,p 2 ) = J - 

T(z ±2 ,Az ±1 ;p 1 ,p 2 ) 

with A = to^i^2^3^4 and 

(5.22) N{t_;p 1 ,p 2 ) = r— r ~ 4 ~ 7—, 

so that 

X f dz 

(5.23) — / A(z;t; Pl ,p 2 ) — = N(t;p h p 2 ) 

2m J z 

is Spiridonov's j22j elliptic Nassrallah-Rahman integral (I1.12J) . valid for \tj\, \pk\ < 1 and 
\pip 2 \ < \A\. 

Fix r G lR<o- We introduce the limiting parameter r G M>o by writing A r (x;r;r) 
and N r (r_; r) for the weight function A(exp(2mrx);t;pi,p 2 ) and the norm N(t;pi,p 2 ) with 
parameters 

tj = exp(27rr(Tj — 1)), pi = exp(27rr/r), p 2 = exp(— 2nr). 
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The conditions \tj\ < 1 and \piP2\ < \A\ then translate to the parameter constraints ()5.19|) 
for the corresponding hyperbolic Nassrallah- Rahman integral (J5.2(Jj) . 
The elliptic beta integral (|5.2H|) can then be rewritten as 

(5.24) / A r (x;r,T)dx = r- l N r (r,r). 

We show in this subsection that the formal limit r \ of ()5.24|) gives the hyperbolic 
Nassrallah-Rahman integral (jl.lOj) . The limit is based on Ruijsenaars' [19, Prop. III. 12] 
observation that the hyperbolic gamma function is a limit case of the elliptic gamma 
function (see ()6.18|) and ()6.19|) for the explicit limit transition in the present notations). 

We first rewrite the weight function A r (x;r;r) and the norm A^ r (r;r) in terms of the 
renormalized elliptic gamma function T r (z; r) (see ()6.18j) ). This yields the expressions 

r r (±2x, hi — ia ± x; rj v/r 

N r (r;r) 



(exp(— 27rr); exp(— 2^))^ (exp(27rr/r); exp(27rr/V)) c 

n <fc< m <4 f r {2i - ir k - ir m ; r) / 5tt 
x — ^— ; — == exp 1 — t 

with a = r + T\ + r 2 + r 3 + r 4 , where we use the notations 

m 

r r (ai, • • • , a m ; r) = JJ IV (a,-; r), 
i=i 

m 

r r (ai ± zi, . . . , a m ± z m ; r) = Jj r r (aj + z,-, aj - Zj\ r) 

i=i 

for products of renormalized elliptic gamma functions. Consequently, ()5.24j) can be rewrit- 
ten as 

, KO _v X[)=o^r{i-iTj±x\T) Uo<k<m<4^r{2i-ir k -ir m ;r) 

(5.25) / -7=^- dx = c r (r) — ~ 4 ~~ , 

J-i T r (±2x, 5? — ia ± x; r) ]^[ - =0 T r (4« — ia + %Tj\ r) 

with 

c r[ r) = 2«P(^-D) 

r (exp(— 27rr); exp(— 27rr)) oo (exp(27rr/r); exp(27rr/r)J 

The limit r \ of the constant c r (r) can be computed using the modularity of the 
Dedekind eta function. 

Lemma 5.6. 

2 

lim Cr(r) = —j=. 
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Proof. The Dedekind eta function rj(a) for a G EI is defined by 

(5.26) rj(a) = (exp(27ricr); exp(27ricx)) oo exp(7rzcr/12), 
and satisfies 

(5.27) T](~) =r ] (a)V^io z . 

Applying ()5.27|) with a = ir and with a = —ir/r, we can rewrite c r (r) as 

2exp(^(l-r)) 



c r (r) 



-t (exp(-2vr/ r) ; exp(-2vr/ r)) (exp(27rr/ r) ; exp(2vrr/r)) 



Since r G M<o we can take the limit r \ in the latter expression, which yields the desired 
result. □ 

We now further investigate the limit r \ of ()5.25|) using ()6.19|) . We associate with 
r the two deformation parameters q = exp(27rzr) and q = exp(— 2m /t), and we use the 
standard convention that q u = exp(27nwr) and q u = exp(— 2niu/r) for u G C. For the 
weight function and norm in ()5.25|) we obtain by ()6.19|) 

nLo f r (i - iTj ± x; r) [1 ± 2ix, a-i±ix; r]^ 



lim = a 



( r ^ ) ^° f r (±2x, 5i-ia± x- r) IlJ=ofa ± r ]°° 

Um llo<fc< m <4 ?r (2j - ir fc - ir m ; r) = ^ U 4 j=0 [ a ~ T i ~ 3 > r ]°° 



rX0 nt=o r »-( 4 ^-^ a + zr i; r ) n <fc<m<4[ r fc + r « - !;tI°o 



with the constants a and (3 given by 

a = g-|gl g -H(^^) 2 + (^ + |- a ) 2 -^o(^ + |- J ) 2 } ; 

/3 = g-^g^gi{^0<fc<m<4(27 + |- T -fe- T ™) 2 -E, 4 =o(^ + |- a + T - J ) 2 }_ 

Writing out the squares leads to 

li / TTt \ 

(5.29) (3 = q~^q ^ exp ( — j a. 

Combining Lemma f5.fi| ()5.28j) and ()5.29j) . we can formally take the limit r \ in ()5.25j) . 
leading to 

^ 3Q ^ r [1 ± 2ix, a-A±tx; r]^ ^ _ 2q-^q h nj= [ a " r J " 3 ? r ]°c 



nJ=o[ T J ± ^ T ]oo exp(-7rV4) V /3 ^ rio<fc<m<4[ r fe + r m - 1; t] 

Changing the integration variable and using 

exp(— 37ri/4)\/— t = —\J—ir 

we arrive at 



[l±2x,a-4±x;r] oc 2q-2iq2i U j= oi a ~ T i ~ 3 ! r 
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which is the hyperbolic Nassrallah- Rahman integral (jl.lOj) . 

6. Appendix: The hyperbolic gamma function 

In this section we discuss the r-shifted factorial [z; r}^ (see f)1.7p ) and its connection with 
Ruijsenaars' |19] hyperbolic gamma function. For detailed proofs we refer to Ruijsenaars' 
papers [TH] and [2U] . 

We start by recalling the definition of Ruijsenaars' hyperbolic gamma function. We write 

C± = {z G C | Re(z) ^ 0} 
for the open right/left half plane. The integral 

dy ( sm(2yz) z 



(6.1) g(z) = g(a + ,a_;z) . 

J y \2smh(a+y) smh(a_yj a + a^y 

defines an analytic function for (a+, a_, z) G V, where 

V = {{a + ,a_,z) G C x3 |Re(a±) > 0, \lm(z)\ < ^(Re(a+) + Re(a_))}. 

Ruijsenaars' 1 1 9 J hyperbolic gamma function is now defined as the analytic, zero-free func- 
tion on T> defined by 

(6.2) T h (z) = T h (a + , a_; z) = exp(ig(a + , a_; z)). 
The hyperbolic gamma function satisfies the functional equations 

/n q\ T h {z + ia±/2) 

(6.3) — ; — j- = 2 cosh{nz/a T ) 

i h {z — ia±/z) 

whenever the left hand side is defined. Hence I\ extends to a meromorphic function in the 
domain ( ; Qj— j Z ) G C* 2 x C, which we again denote by T h (a + , a_; z). 

The functional equations (|6.3|) imply that the zeros of z \— > T h (a + , a_; z) are located at 



and the poles are located at 



Z >o + \ J ia + + ( z >o + ^ I id- 



1>>o + - J ia + - ( Z> + ~ I ia_. 



2/ V 2 

The zeros and poles are simple when a + /a_ G" Q>o- 
The explicit integral expression for g implies 

(6.4) r\(a + ,a_;z) = T h (a_,a + ;z 



(6.5) T h (ra + ,ra_;rz) = T h (a + ,a_;z), r G JR >0 , 
as well as the reflection equation 

(6.6) T ft (a+, a_; ,z)I\(a + , a_; -z) = I. 
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Note that by the translation invariance ()6.5|) . the hyperbolic gamma function I\ essentially 
only depends on the quotient a + /a_ of the two deformation parameters a±. 
Some special values of I\ are easily computed. Obviously, we have 

(6.7) F h (a+, a_; 0) = 1. 

Applying the two functional equations and the reflection equation, we furthermore have 

_, / ia + m_\ / ia+ ia~\ sinh(7rz/a_i_) 
V h [a + ,a-;z — I — — Fh a+, a_; —z — + 



2 27 V 2 2 J sinh(7rz/a_)' 

Taking the limit z we conclude that 



(6.8) 



~2 2 




(to see that the branch of the square root is the right one, note that r^(a+, a_; x) > for 
a± G 1R >0 and x G iR in view of (jgHJ) and flO} )- 

The hyperbolic gamma function can be explicitly expressed as quotient of trigonometric 
gamma functions when Im(a + /a_) ^ 0. This expression was first obtained by Shintani 
[22], see also [201 Appendix A]. 

Proposition 6.1. Let a± G C + with Im{a + /aJ) > 0. Then 

(— exp(7ua + /a_) exp(— 2nz/a_); exp(27ria + /a_)) oo 



T h (a+,a_; 2) 



(— exp(— 7ria_/a + ) exp(— 2nz/a + ); exp(— 27ria_/a + )) ( 
ni f a + cl-\\ ( iriz 2 



x exp — — 1 exp 



24 V a_ cl+ J J V 2a + a 



Proof. We sketch a proof because the infinite product expression for plays such a crucial 
role in the present paper. Write T h (z) for the right hand side of the desired identity. It is 
easily verified that Th(z) is meromorphic in z, having the same poles and zeros as Th{z). A 
direct check shows that Th(z) satisfies the same functional equations as r^(^). Thus 
Th/Fh is an entire, bounded function, hence a constant. The constant is one since 



ia_ za 4 
"2 2 




which follows from the modularity (|5.27jl of the Dedekind eta function ()5.26|) . □ 

Note that an infinite product expression for T/ l ( ) with Im(a+/a_) < can be 

directly derived from Proposition Id II bv applying (jfi.4j) . 

Using Proposition I6.ll the rather harmless looking reflection equation ()6.6|) turns into 
a nontrivial infinite product identity. This identity can be proven without referring to 
the integral representation of F^ using Jacobi's triple product identity, Jacobi's inversion 
formula and the modularity of the Dedekind eta function. 
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Ruijsenaars [20J established rather delicate asymptotic bounds for the hyperbolic gamma 
function. We formulate here a weaker version of these bounds which suffices for our pur- 
poses. It can be stated as 

T h (a + ,a-;z) = O ^exp i^r ^a )) ' ±0 °' 

uniformly for lm(z) in compacta of K. and for a± in compacta of C+. The precise meaning 
is as follows, cf. fK)\ Thm. A.l]. 

Proposition 6.2. Let K± C C + and if C R k compact subsets. There exist positive 
constants R = R(K + , K_; K) and C = C(K + , K_; K), both depending on K± and K only, 
such that 



T h (a + ,a_;z) exp 
when Im(z) £ K and a± G K± . 



ITTZ 2 



2a+a_ 



< C, Re(z) ^ R 



The precise connection between and the r-shifted factorial [^t^ (see f)1.7|) ) is as 
follows. Fix r G C_ fl H and write q = exp(27rir) and q = exp(— 2ni/r). Recall the 
notational convention q u = exp(27rzru) and q u = exp(— 27uu/r) for u G C. The r-shifted 
factorial [z; r]^ (see (|1.7|) ) can then be expressed in terms of the hyperbolic gamma function 
by 

(6.9) [z- r]oo = qi*- 1 *-^) ^q~Uq^r h (~, 1; i (z - \ - j . 

Indeed, substituting the infinite product expression for Th from Proposition 16. II in the right 
hand side of ()6.9j) gives precisely the expression (jl.7j) defining the r-shifted factorial. Work- 
ing with the r-factorial [z; r]^ has the advantage that formulas have a similar appearance 
as in the usual basic hypergeometric (trigonometric) setup. The disadvantage is the loss 
of symmetry in the parameters a± (see (I6.4|) ). We end this section by reformulating the 
above properties of the hyperbolic gamma function in terms of the r-shifted factorial. 

Formula ()6.9|) shows that [z;t]oo has a meromorphic continuation to (z,r) G C x C_, 
which we again denote by [z; r}^ (in fact, in view of (jl.7|) and (|6.4j) . [z; r]^ extends to a 
meromorphic function in (z, r) G C x (C \ R>o)). The functional equations (|6.3|) become 

[2 + 1;t]oo = (l-^^rjoo, 

(6.10) . 

[z-t ;r] oc =(l-g exp(-2mz)) [z; t]^. 

The zeros of z t— > [z; r]^ are located at 

(6.11) -Z <0 + Z >0 , 

T 

and the poles are located at 

(6.12) -Z >0 + Z< . 

T 
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The zeros and poles are simple when r (jL Q<o- Property ()6.4j) translates into 
(6.13) [z;r- 1 ] 00 = [-- + - + l-r] 00 

T T 

for t G C_. The reflection equation ()6.6|) becomes 

(6-14) [- + — + z, - + — - z; rjoo = g 24g 24g 2 

for r G C_. Finally, the asymptotic bounds for the r-shifted factorial, deduced from 
Proposition 16. 2\ become 



(6.15) [z- tU = O {q\ z -*--r) / 2 J , Im(z) -> oo, 

uniformly for Re(z) in compacta of R and for r in compacta of C_, and 

(6.16) [z; rjoo = 0{1), Im(z) -> -oo, 

uniformly for Ke(z) in compacta of R and for r in compacta of C_. Using the more precise 
asymptotic estimates for the hyperbolic gamma function in [19, Thm. A.l], we actually 
have the limit 

(6.17) lim [z;t] 00 = 1 

Im(2)-»-oo 

uniformly for Ke(z) in compacta of R and for r in compacta of C_. 

We end this section by rewriting the hyperbolic degeneration (see Prop. III. 12]) of 
Ruijsenaars' elliptic gamma function in our present notations. Ruijsenaars' elliptic gamma 
function G(r,a + ,a_; z) (see e.g. |T[JJ Prop. III. 11]) relates to the elliptic gamma function 
(dll) by 

G(r, a + , a_; z) = r(exp(2ir^ — a+r — a_r); exp(— 2a + r), exp(— 2a_r)) 

for Re(a±r) > 0. Thus 23 Prop. III. 12] becomes 

_/ / N \ { 2irz — nia^ — ixia_ 

lim Tyexp(27rirz); exp(— 27ra + r), exp(— 27ra_r)J exp 



r\o ' \ 12ira+a- 

ia+ %a_ 



for a± G R>o- We now take a + = — - and a_ = 1 with r G R<o, and we denote for 
simplicity 

(6.18) T r [z;r) = r(exp(27rir,2); exp(27rr/r), exp(— 2nr)) exp f 

By ()6.9|) . we arrive at the limit 

limf r (2:;T) = q~(^-^ lz ) /4 g3sg [r -1 - iz; 

r\0 
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with q = exp(2nir) and q = exp(—2iii/r). Applying the reflection equation (j6.14|) leads 
to the limit 

(6.19) limf r (^;r) = / 4 g -sgs[i + iz;r]^ 

r\0 

for r G R <0 . 
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